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Abstract 

We derive possible corrections to the statistical predictions of quantum mechanics in measure- 
ment over ensemble of identically prepared system based on a hidden variable model of quantization 
developed in the previous work. The corrections are characterized by a dimensionless parameter 
a and the prediction of quantum mechanics is reproduced in the formal limit a — > 0. Quantum 
mechanics is argued to be reliable for sufficiently low quantum number. 
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I. INTRODUCTION 



For almost nine decades since its completion, quantum mechanics has been claimed to 
be confirmed by a wealth of experimental tests with unparalleled accuracy. Nevertheless, 
in view of its operational character the yet unsettled foundational problems j2|, and 
the difficulties in unifying quantum mechanics with general relativity 3], it is imperative 
to ask whether quantum mechanics is an accurate approximation of a deeper theory {4- 
Q. There are at least two possible directions to exercise this question. First, one keeps 
the formalism of quantum mechanics unchanged as far as possible and speculates a small 
correction to the fundamental equation of the theory. The other direction is to construct a 
conceptually new theory which reproduces the empirical statistical prediction of quantum 
mechanics as certain limiting case. In this latter approach, the operational formalism of 
quantum mechanics should be shown as emergent. Obviously, this approach is preferable 
to attack the foundational problems of quantum mechanics and the problems of quantum 
gravity. 



On the other hand, in our previous work 



17] , we have proposed a hidden variable model 



for quantization by modifying the classical dynamics of ensemble of trajectories parame- 
terized by a hidden random variable. We showed for a wide class of important dynamical 
systems that, given the classical Hamiltonian, the modified equations can be put into the 
Schrodinger equation with a unique Hermitian quantum Hamiltonian and Born's statistical 
interpretation of wave function, if the probability density (mass) function of the hidden 
variable A takes a specific form as 

P Q (\ ] h) = ±5(\-h) + ±5(\ + h). (1) 

Namely, A is an unbiased binary random variable which can only take values ±h. We 
have also shown that in this case one can always identify an "effective" velocity field which 
numerically is equal to the "actual" velocity field of particle in pilot-wave theory [3]. This 
then allows us to follow the description of pilot-wave theory on measurement without wave 
function collapse. Yet unlike pilot-wave theory, our model is inherently stochastic and the 
wave function is not physically real field. 



In the present paper, we shall further elaborate the hidden variable model of Ref. 



17| 



in the case when the distribution of A is allowed to deviate from Eq. (CQ) while satisfies the 
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following general condition: 



P(X;h,a) = P(-X;h,a), 
with A ^0, and lim P(X; H, a) = P Q {\; K). (2) 

Hence, A is non-vanishing (P(0; h, a) = 0), its distribution function is even so that it is 
unbiased, and we introduce a new parameter a so that Eq. ([1]) is recovered as a formal limit 
o — > 0. This will then be shown to offer possible corrections to the prediction of quantum 
mechanics in measurement over ensemble of identically prepared system. The discussion will 
be facilitated by considering concrete models of measurement of angular momentum. First, 
we shall show that there is an inherent broadening of spectral line which is purely induced 
by the distribution of the hidden variable. Accordingly, the Born's statistical rule will also 
be shown to have small correction characterized by a. 



II. HIDDEN VARIABLE MODEL FOR QUANTIZATION: QUANTUM HAMIL- 
TONIAN FOR ANGULAR MOMENTUM MEASUREMENT 

Let us consider the dynamics of iV particles with configuration coordinate q = 
(qi, q 2 , . . . , ^at). The classical Hamiltonian is denoted by H_(q,p; t), where {p.} is momentum 
canonically conjugate to {q{\ and t is time. All mathematical symbols with "underline" will 
be used to denote physical quantities satisfying the law of classical mechanics. The classical 
dynamics of the particles then follows the Hamilton- Jacob i equation 

d t S(q;t)+H(q,d g S(q;t);t) = 0, (3) 



where S_(q;t) is the Hamilton principle function (HPF) so that p = d q S_ [19|]. To solve the 
above equation, one needs to set up an initial HPF S_(q; 0) which implies an initial classical 
momentum field p(q; 0) = d q S_(q;0). A single trajectory in configuration space is picked up 
if one also fixes the initial position of the particles. 

Let us then consider an ensemble of classical system so that the probability density 
function of the position of the particles is denoted by p{q;t). It must satisfy the following 
continuity equation: 

dtp + d q ■ (v(S)p) = 0, (4) 
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where v — . . . ,v N ) is the classical velocity field. In the above equation we have made 
explicit the possible dependence of the classical velocity field y_ on the HPF S_, which can be 
obtained from the classical Hamiltonian through the Hamilton equation: 



f = MS), (5) 

iP=d qi S} 



dH 

dp . 
—i 

where i — 1, . . . , N are some functions determined by the choice of H_ as above [20]. The 
dynamics and statistics of the ensemble of classical trajectories are then obtained by solving 
Eqs. (j3J), dH) and (jSJ) in term of S(q; t), p(q; t) and v(q; t). 

Now let us develop a general scheme to modify the above classical dynamics of ensemble 
of trajectories 1^]. To do this let us introduce two real- valued functions S(q, A; t, H, a) and 
Q(q,X;t,h,a), where A is a hidden random variable whose probability density function is 
assumed to satisfy Eq. (J2J). They are supposed to take over the role of S_ and p in the 
modified dynamics. Hence, Q(q, A; t, h, a) is the joint-probability density that the particles 
are at configuration space q and the value of hidden random variable is A. The marginal 
probability densities of the fluctuations of q and A are thus given, respectively, by 

p(q;t, h,a) = J d\Q(q, \;t,h,a), 
P(\-H,(t) = J dqtt(q,\;t,h,a). (6) 

In the following, for notational simplicity, we shall sometime not make explicit notationally 
the dependence of any functions on h and a. 

Next, let us proceed to assume the following general rule of replacement to modify Eqs. 
© and © [l7|: 

p i — y f2, 

d qi S^d gi S+^^, z = l,...,iV, 
dtS^d t S+^ + ±d q .f(S), (7) 

where the vector-valued function / = (/i,...,/at) is defined by Eq. (jSJ). Let us first 
show that the replacement of Eq. (JTj) possesses a consistent classical correspondence if 
S — > S_ so that the Hamilton- Jacobi equation of d2J) is restored (notice that we have used 
the symbol to denote replacement and "— >" to denote a limit). First, using the last 
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two equations of (J7|), for sufficiently small At and Aq = {Aqi}, then expanding AF = 
F(q + Aq; t + At) — F(q; t) m d t FAt + d q F ■ Aq, for any function F, one has 

AS^AS+±(^ + d q .f(S)At). (8) 

One can see that in the limit S — > S_, in order to be consistent then the second term on the 
right hand side has to be vanishing. Namely one has dVt/dt = —Qd q ■ v_, by Eq. (jSJ). This is 
just the continuity equation of (HI). In other words, in the limit S — > S, p(q; t) = J d\fl has 
to approach p, p — > p. 

Let us apply the above modification of classical mechanics to a class of von Neumann 
model of measurement of angular momentum. A different model of measurement will be 
given in Section III I CI To do this, let us consider the dynamics of two interacting particles 
with coordinate q = (91,92)- The first particle represents the system whose angular mo- 
mentum is being measured and the second particle represents the measuring apparatus. To 
make explicit the three dimensional nature of the problem, let us put 91 = (x\,y\,Z\). For 
simplicity let us consider the measurement of z— part angular momentum of the first particle 

where p is the conjugate momentum of x\ and so on. 

Next, let us choose the following measurement-interaction classical Hamiltonian: 

Ki = gL zl p 2 = g{x\v yx - vxp^Ptf ( 10 ) 

where g is the coupling parameter. Let us assume that the interaction is impulsive so 
that one can ignore the free Hamiltonian of each particles. First, L Zi is time-invariant: 
dL /dt = {L , Hf} = where {•, •} is Poisson bracket. On the other hand, one also has 
dq^jdi = {q2,ILi} = 9L. Zl - Integrating, one thus obtains 

q 2 (t) = q 2 (0) + gL Zi t. (11) 

Hence, one can infer the value of L_ of the system (first particle) prior to the measurement 
from the initial and final positions of the apparatus (second particle). 

Let us now consider an ensemble of identically prepared angular momentum measurement 
and investigate the modification imposed by Eq. (J7|) to the classical dynamical equations 
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that govern the ensemble of trajectories. To do this, first, given the classical Hamiltonian 
of Eq. fflUj) . the Hamilton- Jacobi equation of ()3]) becomes 



d t S + g(x 1 d yi S - yi d Xl S)d q2 S = 0. (12) 

On the other hand, substituting Eq. (fTOj) into Eq. fl5]), the classical velocity field is given by 

v Xl = -gyid^S, v yi = gxid q2 S_, v Zi = 0, 

v 2 = g(x 1 d yi S- yid xl S). (13) 

The continuity equation of (@J then becomes 

dtp - gyid Xl (pd q2 S) + gx x d yi (pd q2 S) 

+gx 1 d q2 (pd yi S) - gyid q2 (pd Xl S) = 0. (14) 

Next, from Eq. (1131) and the definition of / given by Eq. (jSJ), one has 

f xl (S) = -gyid q2 S_, f yi = gxid q2 S_, f Zl (S) = 0, 

f 2 (S) = g(x 1 d yi S-y 1 d Xl S), (15) 

so that d q ■ f(S) = 2g(xid q2 d yi S — y\d q2 d Xl S). Substituting this into Eq. (J7J), one then 
obtains 

p i — y ^2, 

A d xl Vl 



d xl S^d xl S + 



2 n ' 



d t S ^d t S + --^- + gX( Xl d yi d q2 S - Vl d Xl d q2 S). (16) 

Let us proceed to see how the above set of equations modify Eqs. (ff2l) and (fl4l) . Imposing 
the first four equations of ( {TBI) into Eq. ( fl4|) one obtains, after a simple calculation 

d t tt - gyid Xl (Qd q2 S) + gxid yi (ttd q2 S) + gxxd^iVtd^S) 
-gyidq 2 (Md Xl S) - g\(yid Xl d q2 Vt - xid yi d q2 ft) = 0. 

(17) 
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On the other hand, imposing the last four equations of (Tl6|) into Eq. (fl2|) . one has, after an 
arrangement 

d t S + g(x 1 d yi S - yi d Xl S)d q2 S - g X 2 {x x ^^ 

-i/i ^|^) + ±(m - gyi d xl (Qd q2 s) 

+gx 1 d yi (Qd q2 S) + gxxd q2 (Qd yi S) - gyid q2 (Qd Xl S) 

-g\(yid Xl d q2 tt - xid yi d q2 tt)^) = 0, (18) 



where we have defined a real- valued function R = \/Q and used the following identity 
(l/4)(^^.fi/n 2 ) = (l/2)(d qi d qj n/n)-(d qi d gj R/R). Substituting Eq. (JEJ) into Eq. (QSJ), 
the last term in the bracket vanishes to give 

d t S + g(x 1 d yi S - y 1 d Xl S)d q2 S 

-gX\ Xl d -^- yi d -^)=0. (19) 

The dynamics of ensemble of trajectories is then determined by pair of coupled Eqs. ffT7|) 
and (Tl9l which depend on the random hidden variable A. 
Now let us assume that Q satisfies the following condition: 

n(q,X;t) = n(q,-X;t), (20) 

so that -P(A) = f dqQ(q, \;t) = P{— A) as required by Eq. (j2J). In this case, S(q,X;t) and 
S(q, —A; t) satisfy the same differential equation of ffl9|) . namely the last term on the left hand 
side is insensitive to the sign of A. Hence, assuming that initially S(q, A; 0) = S(q, — A; 0), 
one obtains 

S(q,X;t) = S(q,-X;t). (21) 

This can be used to eliminate the last term on the left hand side of Eq. ( fT71) . That is, taking 
the case when A is positive add to it the case when A is negative and divided by two one 
gets 

d t n - gy!d Xl (VLd q2 S) + gxid yi {Vtd q2 S) + gxid q2 (ttd yi S) 

-gyid q2 {Vtd Xl S) = 0. 

(22) 
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Recalling that A 7^ as required by Eq. (J2J), let us further define the following complex- 
valued (wave) function: 

A;t) = VOexp (z— -J = i2exp (i-^ -A, (23) 

where we have notationally omit the dependence of W on H and a. From Eq. (|6]), the 
probability density for the position of the particles is thus 

p(q;t)= [ d\\V\ 2 . (24) 



Equations (TH?]) and (I22p can then be rewritten into the following generalized Schrodinger 
equation: 

A 2 

i\\\d t V = -g\ 2 (x 1 d yi - yid Xl )d q ^ = g—L Zl p 2 q, (25) 

where L Zl = —iH(xxd yi — y\d xi ) and p 2 = —ihd q2 are the quantum mechanical z— angular 
momentum and linear momentum operators pertaining to the wave functions of the first 
and second particle, respectively, and we have assumed that the spatiotemporal fluctuations 
of A is ignorable as compared to that of S. 

Let us consider a specific case where Q is separable Q(q, A; t, h,a) = p(q; t, h, cr)P(A; h, a) 
and the distribution of A is given by Eq. as assumed in Ref. [17| . Further, let us define 
a new complex- valued function 

*q(?;*) = VpS*)^^), (26) 
where Sg(q; t) = S(q, ±h; t). Then, Eqs. ( )25|) reduces into the Schrodinger equation 

im*Q = Wq, (27) 

with quantum Hamiltonian Hi 

Hi = gL Zl p 2 . (28) 

Equations f l27|) and f )28|) are the model employed by von Neumann to discuss quantum mea- 
surement 21f| . The above result can be extended to the measurement of angular momentum 
along the x— and y— directions by cyclic permutation of (x, y, z). In this case, L Zl in Eq. 
( 12 8 p is replaced by L Xl and L yi , the quantum mechanical angular momentum operators along 
the x— and y— directions, respectively. We have thus reproduced the results of canonical 
quantization as a specific case of our hidden variable model. We have also shown in Ref. [l^] 
that unlike canonical quantization, the above method of quantization is free from operator 
ordering ambiguity. 
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III. POSSIBLE CORRECTIONS TO QUANTUM MECHANICAL PREDICTIONS 



In this section we shall go beyond quantum mechanics by assuming that the distribution 
of the hidden variable A satisfies Eq. (j2j) rather than Eq. (pfj) . We thus have to start from the 
generalized Schrodinger equation of (I25|) . Various possible corrections to the prediction of 
quantum mechanics will be given. In general, the prediction of standard quantum mechanics 
will be argued to be reliable only for sufficiently low quantum number. 

A. Hidden random variable induced broadening of spectral line 

Let us discuss measurement of angular momentum in ensemble of identically prepared 
system so that the initial wave function of the system (first particle) ip{q\) is given by one 
of the eigenfunction of the angular momentum operator ip(qi) = (j>i(qi), L Zl (j>i = l(f>i, where I 
is the eigenvalue. Further, let us denote the initial wave function of the apparatus (second 
particle) by v^ofe), assumed to be sufficiently localized. The total initial wave function of 
the system-apparatus is thus given by 



We have thus made an idealization that the initial wave function is independent of A. Recall 
that in this case, according to the standard quantum mechanics, each single measurement 
event will give outcome / with certainty (probability one). This is one of the postulate of 
quantum mechanics. 

Let us solve Eq. (1251) with the initial condition given by Eq. (1291) . To do this, let us 
assume that after interval time-span t of measurement-interaction, the wave function can be 
written as 




(29) 



V(q,\;t) = 4> l (q 1 )<p(q2,\;t). 



(30) 



Inserting this into Eq. (l25l) and keeping in mind that L Zl <pi = Ify, one has 



d t (p + gl'd q2 (p = 0, 



(31) 



where I' 



depends on A as 




(32) 
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Equation (IHT1) can then be directly integrated with the initial condition <p(q2, A; 0) = 
(p (q 2 ) to give 

(p(q 2 ,X;t) = (f (q 2 - gl't). (33) 
Inserting this back into Eq. fj30|) . one has 

A; t) = MQi)MQ2 - g\\\lt/h). (34) 

Hence, in each single measurement event, the wave function of the apparatus becomes cor- 
related to the initial state of the system and is shifted an amount of gl'(X)t. This means 
that at the end of each single measurement event, the initial position of the second particle 
(the apparatus pointer) is shifted uniformly as 

q 2 (t,\) = to(p)+gl'(\)t. (35) 

Now let us interpret the above formalism in similar way as with classical measurement. 
As discussed in the previous section, in the latter case, after time-span of measurement- 
interaction t, the position of the apparatus-particle is shifted as q 2 {t) = q 2 (0) + gL Zi t. From 
this, one infers the result of measurement to be given by L . Similarly, it is natural to 
interpret Eq. (I3"5"j) that the outcome of each single measurement event is given by /'(A) = 
\X\l/H of Eq. (|32|) . Here we have applied the result shown in Ref. [17j that it is possible 
;o probe the pre-existing value of the initial and final positions of the apparatus particle 



221 ] . Hence, instead of obtaining a sharp value / as postulated by the standard quantum 
mechanics, one obtains a random value /'(A) which depends on the value of the hidden 
variable A. One can also see that when the distribution of A is given by Eq. ([TJ so that 
A = ±h, then the randomness of the outcome of single measurement disappears and one 
regains the prediction of quantum mechanics: l'(±h) = I with probability one. For general 
distribution of A satisfying Eq. fl2]), we have thus a random correction to the prediction of 
quantum mechanics: even when the initial wave function of the system is given by one of 
the eigenfunction of the angular momentum operator, the result of each single measurement 
will still be random with statistical properties determined by the distribution of A. 
Hence, given the value of I, the probability density to get I' is 

h 



P{l'\l) = -(f + (\;h,a)\ x ^ h + f4\;H,a) 



2jf + (X;h,a)\ x= ^ (36) 
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where f + (X;H,a) and f-(X;h,a) are part of P(X;h,a) defined on the positive and nega- 
tive axis of A respectively, and in the second equality we have used the assumption that 
f + (\;h,cr) = f_(—X;h,a) of Eq. ([2]). In the limit o — )■ one has lim^o /+(A; h, a) = 
(1/2)5(X — h) by Eq. ([T]) so that we reproduce the prediction of quantum mechanics 

KmP(l'\l) = ^6(j(l>-l))=5(l'-l), (37) 

that is, in each single measurement event, one always obtains /' = /, as expected. 

Let us proceed to discuss the statistical properties of /' in term of the statistical properties 
of the hidden variable A, P(A; h,a). Recall that P(A; h,a) must satisfy Eq. (j2J). There are 
then infinitely many P(A; H, a) fulfilling this requirement. Let us give a general method 
to construct such probability density function. First, since P(A; H, a) = P(— A; H, a) then 
it is sufficient to fix the form of P(A; h, a) on the half line A > 0. Let us then pick up a 
non-negative function denoted by P + (X;h,a) which is defined on A > 0. Further, let us 
assume that P+(A; h, a) is normalizable, J °° dXP + (X; h, a) = 1, and possessing the following 
limiting property: 

UmP + (X;h,a) =5(X-h). (38) 

a thus measures the width of P + (X;h, a). The desired probability density for the hidden 
random variable A can then be constructed as 

P(A; h, a) = l -P + {X; h, a)U(X) + \p+{~\ K, v)U(-X), (39) 

where U(X) is the Heaviside step-function, namely U(X) = 1 for A > and U(X) = for 
A < 0. It is then evident that the so-constructed P(A; H, a) possesses the required symmetry 
property P(A; h, a) = P(— A; h, a). Moreover, in the formal limit a — > 0, one obtains, by the 
virtue of Eq. flUHJ 

lim P(A; h, a) = h(X - H) + ^(A + h) = P Q (X; H), (40) 

as required by Eq. (J2J). The prediction of quantum mechanics is thus regained in the limit 
of vanishing a which is equal to the vanishing of the width of P+(A; h, a). 

The mean and variance of the fluctuation of I' conditioned on the value of I (quantum 
number) can then be expressed as follows. First, given I, the mean of I' is 

M 1 [P(l'\l)] = -J dX\X\P(X;h,a) = -M x [P + {X-h,a)). (41) 
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Notice that in the limit o — > 0, one has Mi[P(A; ft, a — > 0)] = ft by Eq. ( )38|) . so that one 
regains the prediction of quantum mechanics Mi[P(l'\l)] — > I. In general, for non- vanishing 
a, however, Mi[P + (X; ft, a)} ^ ft and is independent of I so that there is a correction to the 
prediction of quantum mechanics which is proportional to I (the value predicted by quantum 
mechanics). Similarly, the second moment is given by 

M 2 [P(l'\l)] = jpJ dXX 2 P + (X ] h,a). (42) 

The variance of /' given the value of I is thus 

Var[P(/'|/)] = ^Var[P + (A;ft,<7)]. (43) 

Again in the limit a — > 0, one regains the prediction of quantum mechanics Var[P(/'|/)] — > 0, 
by Eq. (138]) . Hence, in general for non- vanishing a, there is a finite broadening of the spectral 
line given by the width of P+(X; h, a) and is proportional to /. 

Let us take a concrete statistical model by assuming the following form of P+(X; h, a): 

It is the log-normal distribution with location parameter In h, scale parameter a and thus 
mode (the position of its maximum) Am = hexp(— a 2 ); x = In A is normally distributed with 



mean In h and width a 23J]. Hence, in the limit of a — > one indeed has lim CT ^o P+{X; h, a) 



5(X — h), as required by Eq. ( |38l) . The mean, second moment and variance are given by 

Mi[P + (A; h, a)] = he' 7 ^ 2 , M 2 [P+(A; ft, a)] = h 2 e 2 ° 2 

Var[P+(A; ft, a)] = ftVV' ~ !)• (45) 

In this case, the profile of the broadening of the quantum mechanical spectral line with 
quantum number I can be obtained by inserting Eq. ( |44"1) into Eq. (1361) . noticing /+(A; ft, cr) = 
P + (A; ft, cr)/2, to give 

where /'// > 0. The profile of the broadening is given by the log-normal function and thus 
not symmetric. Next, the conditional average of V given in Eq. (1411) is 

2 

M 1 [P(r\l)] = le°*l 2 + + 0(a 4 ). (47) 
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Hence, the absolute value of the average of I' is always larger than the prediction of quantum 
mechanics. The correction to the prediction of quantum mechanics is proportional to the 
latter and also to the value of a 2 . Further, the second moment of V conditioned on the value 
of I is given by M 2 [P(l'\l)] = l 2 e 2a2 . The variance of V given the value of I is thus 

Var[P(/'|/)] = ZVV 2 - 1) « ZV + 0{a A ). (48) 

In the quantum limit where o — > we regain the prediction of quantum mechanics: 
M l [P{l'\l)\ I, M 2 {P{l'\l)\ I 2 and Vsx[P{V\l)\ 0. 

The above results suggest that the prediction of quantum mechanics is reliable only for 
sufficiently low quantum number Namely, the deviation from the prediction of quantum 
mechanics grows as ]/] increases. In particular for a symmetric log-normal model with a given 
cr, the prediction of quantum mechanics is ambiguous for |Z| satisfying \l\a 2 /2 « A/, where 
Al is the quantum mechanical spectral spacing. In general statistical model, for sufficie ntly 



large value of |Z|, Var[P(/'|/)] ~ (A/) so that the quantum discreteness is smoothed out [24]. 



B. Modified Born's statistical rule 

Now let us consider the general case when the initial wave function of the system ip(qi) 
is not necessarily the eigenfunction of the angular momentum operator. To do this, first, 
notice that Eq. ff25l) is linear with respect to Hence, since <fii(qi)<Po(q2 — satisfies 
Eq. (|25|) as shown in the previous subsection, their linear superposition over all possible 
values of / 

A; t) = Ci0i(?i)Vo(?2 - gl\Mt/[h), (49) 
i 

also satisfies Eq. ff25l) . Here {q} is a set of complex numbers to be determined as follows. 
Putting t = 0, one obtains 

*(?;0)= ($>M<zi))po(<fc), (50) 
i 

which is separable and independent of A. Hence, the initial wave function of the system 
alone (the first particle) is given by 

V»(9i) = I>&(?i)- ( 51 ) 
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This shows that c\ is the coefficient of expansion of the initial wave function of the system 
ip{qi) in term of the set of orthonormal eigenfunctions of the angular momentum operator 
{fa}, ci = {<f>i\i)). 

One can then see from Eqs. (|49p and f[5Uj) that at the end of the measurement-interaction, 
the wave function of the apparatus-particle separates into a series of packets, each is corre- 
lated to one of the eigenfunction of the angular momentum operator. Namely, for a given 
value of A, the wave packet of the apparatus is shifted an amount of gl'(X)t = g\X\lt/h. 
Hence, if A is fixed and y^ofe) is spatially localized, then for sufficiently large value of g 

<Pi(.Q2, A;t) = (p (q2 -gl\X\t/h), (52) 

does not overlap with each other for different values of /, and each is correlated to a distinct 
eigenfunction of angular momentum operator, <fii{q\). 

Let us now denote the probability density that the apparatus-particle enters the support 
of the wave packet <pi(q2,\;t) as P Vl . Then, the probability density to get the value I' is 
given by 

P(0 = E P ™ p w. (53) 
I 

where P(l'\l) is the probability density to get I' provided that the apparatus-particle is inside 
the support of tpi which is discussed in the previous subsection and is given by Eq. 

It thus remains to calculate P ipr To do this, first, since for sufficiently large value of g, 
!fi(q 2 ,X;t) in Eq. (T49]) does not overlap for different values of I, then the joint-probability 
density that the first particle (system) is at q\ and the second particle (apparatus) is at q 2 
for a given value of A is decomposed into 

Q(q, A; t) = |tf (<z, A; t)\ 2 « ^ |c,| 2 |^(?i)lV/(&. ^ *)f (54) 

i 

Let us note that when | A | is very small one needs a large value of g to separate <pi for different 
values of I otherwise the above decomposition is not valid. Nevertheless, since P(A; h, a) is 
very small in the regime where |A| h, one can argue that its contribution is ignorable. 

From Eq. (1341) . one can see that the joint-probability density that the first particle has 
coordinate q\ and the second particle has coordinate q 2 inside the support of (pi for a fixed 
value of A is given by 

M 2 |<Mgi)| 2 |^fe,A;t)| 2 . (55) 
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The probability density that the second particle is inside the support of the wave packet <pi 
regardless of the position of the first and second particles and the value of A is thus 

P V i = J dq 1 dq 2 d\\ci\ 2 \(pi{q 1 )\ 2 \(pi{q 2 ,\;t)\ 2 = |q| 2 , (56) 

which is just the Born's statistical rule. 

Finally, inserting Eq. (156]) into Eq. ( 153]) . the probability density to get I' can be calculated 

as 

p{i') = Y,nm\ci\ 2 - (57) 
i 

We have thus a modified Born's statistical rule. Since in the limit a — > 0, P(X; h, a) reduces 
into P Q {\; K) given by Eq. © so that one has P(l'\l) 8{V - I) of Eq. (E7J), then Eq. ([57]) 
reduces into 

limP(Z') = V \ci\ 2 5(l' -I), (58) 

a— >0 ' J 

I 

as postulated by quantum mechanics. 

Next, using Eqs. (l5"Tj) and fj4T|) . the average of /' can be calculated to give 

M 1 [P(l')] = J dl'l'P(l') = ^M Mq . (59) 

where Mq = ^ | | 2 is the quantum mechanical prediction for the average value of Hence, 
there is a correction which depends on the value of a. If the initial wave function of the 
system is (pi, one has Mq = I so that one regains Eq. fHTl) . For the case where -P+(A; h, a) is 
log-normal function given by Eq. ( )44|) . then Mi[P + (X; h, a)} = hexp(a 2 /2) so that one has 

M 1 [P(l')] = M Q e° 2 / 2 w M Q (1 + a 2 /2) + 0(a 4 ). (60) 

Further, using Eq. (|42l) . the variance can be expressed as, after an arrangement, 

Var[P((<)] = VarQ^I^iMl + „'}^M, ( 61 ) 

where VarQ = YliU ~ Mq) 2 \ci\ 2 is the variance predicted by quantum mechanics. Again, if 
the initial wave function of the system is (pi, then VarQ = and Mq = I so that one regains 
Eq. fj43|) as expected. Assuming that -P+(A; h,a) takes the form of log-normal function of 
Eq. (HID, one has 

Var[P(0] = VarQe 2 " 2 + M 2 (f 2 - l)e ff2 

w VarQ(l + 2(j 2 ) + Af2 (T 2 + 0( ( T 4 ). (62) 
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From Eqs. (|38|) . ( )59|) and (ISTI) . one can see that the prediction of quantum mechanics is 
regained in the "quantum limit" a — > 0. 

All the above corrections also a ppl y to measurement of other physical quantities. In 
particular, we have shown in Ref. 17) that the measurement of momentum and position 
can also be treated in the same way as the measurement of angular momentum discussed in 
the present paper, so that the corrections to the prediction of quantum mechanics obtained 
above also apply to the measurement of position and momentum. Since in general for non- 
vanishing a one has M2[P+(X; h,a)} > ft 2 , one can then conclude that the variance of the 
measurement results predicted by the present hidden variable model is in general larger than 
the prediction of the quantum mechanics. Hence, for the case of measurement of position 
and momentum of identically prepared system, we will have an inequality which is in general 
stronger than the Heisenberg uncertainty relation. 

Further, it is also evident from the above discussion that there are two sources of ran- 
domness of the outcomes of a single measurement event. The first kind of randomness comes 
from our ignorance of the initial position of the particles. As discussed above, this type of 
randomness determines the probability that the apparatus particle enters the support of the 
wave packet tpi = v?o(<?2 — gl\^\t / ft) which is given by P Vl = |q| 2 . This is the only type of 
source of randomness asserted by pilot-wave theory [181 ] . The second type of randomness 
comes from the fluctuations of A which determines the probability to get I' provided that 
the apparatus particle is inside the support of (pi, P(l'\l). This latter kind of source of 
randomness of single measurement event thus is inherent in the dynamics. 

Let us further discuss the important role played by A. In Ref. 17| . we have shown that in 
the limit o — > our hidden variable model effectively reproduces the mathematical formalism 



of pilot-wave theory 



18) . In most literatures of pilot-wave theory, position is regarded as 



hidden variable. However, it is through the position of the particle that we experience the 
real world. Noticing this fact, Bell then proposed to regard the wave function, which is 
more hidden to us and is assumed to be physically real, as the hidden variable; and call the 
position as the dynamical beable [25[. In this context, A in our statistical model allows us 
to omit the necessity to assume the wave function as physically real. In our model, the wave 
function is only an artificial calculational tool. Its status as hidden variable is replaced by 
A. Moreover, in contrast to pilot-wave theory which is deterministic, A makes the model 
inherently stochastic. It is also the main ingredient of the quantization processes through 
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which we get the Schrodinger equation with a unique physical interpretation. Hence, it is A 
that distinguishes "quantum-ness" from classical statistical mechanics. 



C. Stern-Gerlach experiment 

To give a concrete example and to show the robustness of the statistical model, let us now 
apply the statistical model developed in the previous subsections to investigate the Stern- 
Gerlach experiment of measurement of angular momentum. As will be clear, Stern-Gerlach 
type of measurement is different from the von Neumann model of measurement discussed in 
the previous subsections. 

First, let us discuss the description of ensemble of the Stern-Gerlach experiment in clas- 
sical dynamics. Let us assume that we have a beam of neutral atom whose center of 
mass coordinate is denoted by q 2 = (x 2 , y 2 , z 2 ), each containing an electron with coordi- 
nate qi = (xi,yx, Zx). The interaction between the atom and the magnetic field of the 
Stern-Gerlach apparatus is thus mainly due to angular momentum of the electron. Let us 
assume that this interaction is impulsive so that one can neglect the free Hamiltonian of 
the system during the measurement-interaction. Further, let us assume that the magnetic 
field is non-vanishing only in z— direction B = (0, 0, B z ), B z = B'z2, where B' is some con- 
stant, namely it is a monotonic function along the z— axis [26]. In this case, the classical 
interaction-Hamiltonian can be approximated to be given by 

eB' eB' 
KsG * 2^^- = 2^^ X ^n ~ ™*J> (63) 

where e is charge of electron, m e is its mass, c is the velocity of light. We have also ignored 
the quantum mechanical spin degree of freedom. 
The Hamilton- Jacobi equation of fl3]) thus reads 

d t S + C(z 2 )(xxd yi S - yi d Xl S) = 0, (64) 

where we have defined ((z 2 ) = ^~~^ z 2- On the other hand, from the Hamilton equation of 
(J5j), one obtains the following classical velocity field 

v Xl = -Cizdyu v m = ((z 2 )xx, v Zl = 0, v Z2 = 0. (65) 

The above set of equations give constraints to the dynamics. Using Eq. ( 1651) . the continuity 
equation of (j3J) becomes 

d t p + C(z 2 )(x 1 d yi p-y 1 d Xl p) = 0. (66) 
17 



Hence, the classical dynamics of ensemble of trajectories during interaction with Stern- 
Gerlach magnetic field is given by solving Eqs. ( 164)) and ( 166]) in terms of S_(q,t) and p(q,t). 
Next, from Eq. (1651) . / defined in Eq. (0) is given by 

fxi(S) = -((z 2 )yi, f yi (S) = C(z 2 )x 1} 

/*(£) = 0, f Z2 (S) = } (67) 

so that one has Y^i®qifi(S) = 0- Equation ((Tj) then becomes 

(68) 

Let us investigate how the above set of equations modify Eqs. (IMj) and (166|) . First, since 
Eq. (166!) does not contain S_, then imposing the first equation of (|68|) . one has 

+ C0 2 )(>i^ - yid Xl Q) = 0. (69) 
Further, imposing the last three equations of ( 1681) into Eq. ( 1641) one has 



dtS + CfaKxidnS-ytf^S) 

+ ^ z ^ x ^ n ~ v&ify) = °- ( 7 °) 

Substituting Eq. (169|) . the above equation becomes 

d t S + C(z 2 )( Xl d yi S - Vl d Xl S) = 0. (71) 

We have thus pair of equations (1691) and (1711) which are similar to its classical counterpart 
of Eqs. (IBB) and (IMj). 

Finally, defining complex-valued function as 

# M = Rexp(iS/h), (72) 

where i? = Vs), Eqs. ( )69|) and ( 17 ip can be recast into the following Schrodinger equation: 

ihd t y M = ((z 2 )L Zl V M . (73) 
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Now let us discuss the detail process of the Stern-Gerlach experiment. Let us first assume 
that the initial wave function of the total system prior to entering the Stern-Gerlach magnetic 
system is separable given by 

*m(<?i, (l2\ 0) = <M<?iVo(<72), (74) 

where again <pi is the eigenfunction of quantum mechanical angular momentum operator L Zl 
corresponding to eigenvalue I: L Zl <pi = Ify. Let us assume that after spending time t inside 
the magnetic system of the Stern-Gerlach apparatus, the wave function is still separable 
given by 

*Af(<7i, q.2\ t) = <M<?iM<?2; t), (75) 
where <f(q 2 ] 0) = (p (q 2 ). Inserting Eq. (175]) into the Schrodinger equation of (1731) one obtains 

ihd t ip = ((z 2 )l(p, (76) 

which can then be directly integrated to give 

tp(q 2 ; t) = ipo(q 2 ) exp(-i[iltz 2 /h), (77) 

where fi = eB'/(2m e c). Inserting this back into Eq. f!75|) . the total wave function at the 
exit of the Stern-Gerlach magnetic system at t = T is thus 

* M { qi , q 2 ; T) = Mqi)M<l2)e-^ AlZ \ (78) 

where we have defined A/ = filT. One can thus see that the atomic wave function gets a 
/—dependent phase with a wave vector along the z— direction given by A;. 

Now, after passing through the Stern-Gerlach magnet, let us assume that the atom is 
free. Thus the time evolution afterward is governed by a generalized Schrodinger equation 

i\\\dMq,\;t) = -^-^ 2 2 Vl/(g,A;t), (79) 

where m a is the mass of the atom and we have ignored the electronic free Hamiltonian 
and the irrelevant x 2 and y 2 part of the atomic wave function. The detail derivation of 
Eq. (179]) by imposing Eq. (J7J) to classical dynamics of ensemble with classical Hamiltonian 
£L = p 2 /(2m a ) is given in [17]. Equation f J79|) has to be solved subject to the initial wave 
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function at t — given by Eq. (!78|) : 0) = ^m(Q] T). To do this, let us take a concrete 
model when the initial atomic wave function is Gaussian 



^(^-exp^-^j, (80) 

up to normalization constant, where <jq is the width of the Gaussian. The Schrodinger 
equation of f!79|) can then be solved exactly to give 

( z 2-^-ttQ 2 a i A; |A| , 

\I/(g; t, A) ~ 4>i{qi)e *«tm*o 1 » 122 aman 1 ^ (81) 

where cr t (A) = o"o(l + z|A|t/2m a o"o). 

One can thus see from Eq. (1HTI) that at time £ (measured just after the the particle 
leaving the magnetic field of the Stern-Gerlach apparatus) , the initial atomic wave function 
is uniformly shifted by an amount 

Oi(t; A = = m±±, g M = (82) 

m a n n m a 

We can thus admit that result of single measurement event is random given by /' = \X\l/h 
reproducing the results obtained using the von Neumann model as discussed in the previous 
subsections. Hence, we can proceed as in the previous subsections to derive various correc- 
tions to the prediction of quantum mechanics by assuming a statistical distribution of A. 
The above scheme can also be straightforwardly generalized to arbitrary initial electronic 
wave function ip(qi). Hence, the precision of Stern-Gerlach experiment conducted so far can 
be regarded to give the upper bound for the assumed small yet finite value of a. 



IV. CONCLUSION AND DISCUSSION 



Following Ref. 17| . we have modified the classical dynamics of ensemble of trajectories 
parameterized by hidden random variable to simulate the quantum fluctuations. The hidden 
variable A is assumed to be non-vanishing, its probability density function satisfies a general 
symmetry condition of Eq. (j2j) so that it is unbiased, independent of time, and is character- 
ized by the reduced Planck constant h and a real-valued dimensionless parameter a. The 
statistical model is then applied to discuss the measurement of angular momentum in two 
different models, that of von Neumann and Stern-Gerlach, over an ensemble of identically 
prepared system. 
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We showed that the prediction of quantum mechanics is regained in the formal limit 
a — > 0, so that for non- vanishing a -C 1 it attains small corrections. First, we showed that 
there is a finite broadening on the spectral line which is purely due to the fluctuations of 
the hidden variable. Accordingly, the Born's statistical rule has to be slightly modified. In 
general, the correction to the prediction of quantum mechanics is larger for higher quantum 
number so that quantum mechanics is reliable for sufficiently low quantum number. These 
results thus allow precision tests of quantum mechanics against our hidden variable model. 
Let us mention that there are reports on precision test of quantum mechanics against possible 



nonlinear modification of Schrodinger equation [27]. One of the interesting feature of our 
model, in this respect, is that possible deviations from the prediction of quantum mechanics 
can be accounted for without giving up the linearity of the fundamental equation. 

Since our statistical model reproduces the prediction of quantum mechanics when a — > 0, 



then in this limit it must necessarily violate Bell's inequality [25l . |28j. In other words, in 
the limit of a — > 0, our statistical model is (Bell)-nonlocal. Since in this case we effectively 
obtain the mathematical formalism of pilot-wave theory, then the source of the non-locality 
can also be argued as due to the presence of a new type of potential, the last term of Eq. 
f fl9|) . which is called as quantum potential in pilot- wave theory. However, unlike pilot- wave 
theory, there is no instantaneous interaction between space-like particles for the effective 
velocity in our model is not the actual velocity of the particles but is an average velocity 

17]- 

An interesting question then arises for the general case when a is not vanishing. Since 
in this case our statistical model suggests a small yet finite correction to the prediction 
of quantum mechanics, then one can ask whether the Bell inequality is still violated or 
not. We expect that this issue is intimately related to the argumentation that precision 
iest of quantum mechanics might give non-trivial limitation to the Bell non-locality test 
29I ] . In this context, it is interesting to note that while there have been many experimental 
tests of quantum mechanics against hidden variable theories in view of Bell nonlocality 
and noncontextuality {30], hitherto, to our knowledge, there is no experiments which aims 
to directly test the precision of quantum mechanics against hidden variable theory. This 
might be partly due to very few testable predictions of hidden variable theories that suggest 



correction to the prediction of quantum mechanics 
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It is then interesting to elaborate further extension of the hidden variable model discussed 
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in the present paper. The first alternative is to give up the demand that the distribution 
of A is symmetric P(A; H, a) = P(— A; h, a) for all values of a while keep assuming that the 
symmetry is restored in the limit a — > 0, lim CT ^ P(A; h, a) = Pq(A; h) = S(X — h)/2 + 8(X + 
h)/2. This will again reproduce the prediction of quantum mechanics as the first order 
approximation due to the smallness of a. Since the distribution of A is no more symmetric, 
rather than using the generalized Schrodinger equation of (|25|) . one should start from pair 
of Eqs. ffTTj) and ffT9l) in the case of von Neumann model. The other alternative is to drop 
the assumption that h and a are constant. Namely we allow one or both of them to depend 
on space and time: h = h(q;t) and a = a(q;t). Note that even if they are, assuming that 
the spacetime fluctuations of H and a is negligibly weak compared to the fluctuations of 
ty(q,\;t), then the generalized Schrodinger equation of (|2"5|) or (|79[) are still approximately 
valid. 
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